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Quantitative theory of the Josephson effect in SFIFS junctions (S denotes bulk superconductor, F — 
metaUic ferromagnet, I — insulating barrier) is presented in the dirty limit. Fully self-consistent numerical 
procedure is employed to solve the Usadel equations at arbitrary values of the F-layers thicknesses, magne- 
tizations, and interface parameters. In the case of antiparallel ferromagnets' magnetizations the effect of the 
critical current Ic enhancement by the exchange field H is observed, while in the case of parallel magnetizations 
the junction exhibits the transition to the yr-state. In the limit of thin F layers, we study these peculiarities of 
the critical current analytically and explain them qualitatively; the scenario of the O-vr transition in our case 
differs from those studied before. The effect of switching between and tt states by changing the F-layers' 
mutual orientation is demonstrated. 

PACS: 74.50.+T, 74.80.Dm, 75.30.Et 
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Josephson structures involving ferromagnets as weak 
link material are presently a subject of intensive study. 
The possibiHty of the so-called "vr-state" (character- 
ized by the negative sign of the critical current Ic) 
in SFS Josephson junctions was predicted theoretically 
0, g |, § |, g 0, §. The first experimental ob- 
servation of the crossover from 0- to 7r-state was re- 
ported by Ryazanov et al. and explained in terms 
of temperature-dependent spatial oscillations of induced 
superconducting ordering in the diffusive F layer. 

More recently a number of new phenomena were pre- 
dicted in junctions with more than one magnetically or- 
dered layer. First, the possibility of the critical current 
enhancement by the exchange field in SFIFS Josephson 
junctions with thin F layers and antiparallel magneti- 
zation directions was discussed in the regimes of small 
S layer thicknesses [10 and bulk S electrodes p^ . 
Second, the crossover to the 7r-state was predicted in 
Ref. for the parallel case even in the absence of the 
order parameter oscillations in thin F layers. Still, the 
physical explanation of these effects and accurate cal- 
culation of their magnitude have not been given so far. 
To make such estimates in the model with thin S elec- 
trodes, one must consider KO-1 type solutions and 
take into account spatial variation of superconducting 
state in the SF bilayers; at the same time, in the bulk 
S case an approximate method was used in Ref. |pl| 
beyond its applicability range This problem is of 

rather general nature, since one may expect from the 
previous knowledge (see, e.g., review 0) that the su- 
percurrent in a short weak link is 77-indepcndent. 
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The above intriguing scenario motivated us to attack 
the problem of the Josephson effect in SFIFS junctions 
by self-consistent solution of the Usadel equations for 
arbitrary thicknesses of the F layers, barrier transparen- 
cies and exchange field orientations. Below we show that 
the O-TT transition in the case of parallel H orientation 
or enhancement of Ic by H in the antiparallel case with 
thin F layers occurs when the effective energy shift in 
the ferromagnets (due to the exchange field) becomes 
equal to a local value of effective energy gap induced 
into a F layer. Under this condition a peak in the local 
density of states (DoS) near the SF interfaces is shifted 
to zero energy. In the models with DoS of the BCS type 
this leads to logarithmic divergency of Ic in antiparal- 
lel case at zero temperature, similarly to the well known 
Riedel singularity of ac supercurrent in SIS tunnel junc- 
tions at voltage eV — 2A. We also describe the general 
numerical method to solve the problem self-consistently 
and apply it for quantitative description of the O-tt tran- 
sition and Ic enhancement in SFIFS junctions. 

The model. We consider the structure of SFIFS 
type, where I is an insulating barrier of arbitrary 
strength. We assume that the S layers are bulk and 
that the dirty limit conditions are fulfilled in the S and 
F metals. Although our method is applicable in the 
general situation of different ferromagnets and super- 
conductors, for simplicity below we illustrate our results 
in the case when equivalent S and F materials are used 
on both sides of the structure (although the directions 
of the exchange field in the two F layers may be dif- 
ferent), both F layers have the thickness dp, and the 
two SF interfaces have the same transparency. At the 
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same time, we do not put any limitations on dp and the 
transparency. 

The Usadel functions G, F obey the normaUzation 
condition + F^FZ^ = 1, which allows the following 
parameterization in terms of the new function <I>: 



F,, 



(1) 



The quantity uj — lu + iH corresponds to the general 
case when the exchange field H is present. However, in 
the S layers H — and we have simply uj — u. 

We choose the x axis perpendicular to the plane of 
the interfaces with the origin at the barrier I. The Us- 
adel equations in the S and F layers have the form 



ojGs dx 



ujGf dx 



9 d 



-A, 



0, 



(2) 
(3) 



where Tc is the critical temperature of the superconduc- 
tors, A is the pair potential (which is nonzero only in 
the S layers), uj is the Matsubara frequency, and the 
coherence lengthes ^ are related to the diffusion con- 
stants D as £,s{F) — •\/i?s(F) /27r Tc. The pair potential 
satisfies the self-consistency equations 



A In 



T 



A — Gs^s sgnu! 



0. 
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In the present paper we restrict ourselves to the cases 
of parallel and antiparallel orientations of the exchange 
fields H in the ferromagnets. 

The boundary conditions at the SF interfaces {x = 
Tdp) have the form |0 (see Ref. for details) 



±7s 



to dx 
^fGf d 



*s = 7- 



uj dx 



bJ ox \ U) LU 



(5) 
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with jb^RbA/pfCf, J = Ps£.s/pf£.f, 



where Rb and A are the resistance and the area of the 
SF interfaces; ps(f) is the resistivity of the S (F) layer. 
At the I interface {x — 0) the boundary conditions read 



3i- 

uji dx 



^Fl = 



d 



^fGfi d 

1B,I — — = (JF2 



UJ2 dx 

$F2 



F2, 



Fl 



LOl dx \ UJ2 L0\ 

with 7s,/ = Rb,iA/ pF^F, 



(7) 
(8) 



where the indices 1,2 refer to the left and right hand 
side of the I interface, respectively. 



In the bulk of the S electrodes we assume a uniform 
current-carrying superconducting state 

Ao exp {i[^ip/2 -f 2mvsx\) 



^{x — =Foo) = 



(9) 



1 + 2Dsm^vl/^uj^ -f 

where m is the electron's mass, Vs is the superfluid ve- 
locity, and is the phase difference across the junction. 

The supercurrent density is constant across the sys- 
tem. In the F part it is given by the expression 



J 



2^ ^ 



d_ 

dx 



d_ 

'd^ 



(10) 



while analogous formula for the S part is obtained if 
we substitute u uj. This expression, together with 
the boundary condition (|^) and the symmetry relation 
F{—uj,H) = F{uj,—H), yields the formula for the su- 
percurrent across the I interface: 

^ - ^^^lni[F*,i-H,)FF2{H2)] (11) 
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[the functions F are related to $ via Eq. (^]. 
The limit of small F-layer thickness: 

min(^i?, \/ F>f/2H). Under the condition 75/7 ^ 1 we 
can neglect the suppression of superconductivity in the 
superconductors. We assume further that the trans- 
parency of the barrier I is small, ^b,i 3> max(l,7B), 
and the SF bilayers are decoupled (the exact criterion 
will be given below). In this case we can set Vg ^ and 
expand the solution of Eq. (||) in the F layers up to the 
second order in small spatial gradients. Applying the 
boundary condition (^, we obtain the solution in the 
form similar to that in SN bilayer ||l^, |l^ : 

'^1,2/'^ 



1,_F2 



with 



1 + lBM^l,2/TrTcGs 

IBM = IbAf/^f, Gs^uj/^l 



Aoexp(Ti(^/2), (12) 



A2 



Substituting Eq. (^2|) into the expression for the super- 
current ( pi] ) we obtain I{tp) = Ic siniy9. 

For the parallel orientation of the exchange fields, 
Hi = H2 = H , the critical current is 

27rT ^ 6'^Gl l-a + n^BAigi 



gRb.i 



(1 - a + iljBAigi) 



4,ag2 
(13) 

where Q, — lo/ttTc, 5 = Ao/ttTc, a — {h^BhiY, h = 
H/ttTc, 51 = 2Gs + 7BA/f^, 92 = {Gs + Jbm^?- 
For the antiparallel orientation. Hi = —H2 
the critical current is given by 

27rr j2 ^^Gl 1 



H, 



eR 



BJ 



O>0 



a/(1 - a -I- Vl^BMgiY + 4ag2 

(14) 
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At h = 1/jBM and small the expression under the 
sum in Eq. (^4|) behaves as l/il, thus at low T the crit- 
ical current diverges logarithmically: /i"' cx hi{Tc/T). 
This effect was pointed out earlier in Refs. [Tl| ]. 

The above results become physically transparent in 
the real energy e representation. Making analytical con- 
tinuation in Eqs. ([^), by replacement uj — > —ie, we 
obtain the expression for the DoS per one spin projec- 
tion (spin "up") NF{e) = IleGF{e) in the F layers 



Npie) 



Re 



(15) 



s = £ + 1bm{£ - H)JAl - eVvrT,, 



which demonstrates the energy renormalization due to 
the exchange field. Equation (|l^) yields Nf{0) = 
R.e(7_BAf/i/\/(7BAf^)^ — 1), which shows that a,t h — 
^llBM the singularity in the DoS is shifted to the Fermi 
level. Exactly at this value of h the maximum of li"^ IS 
achieved due to overlap of two e~^/^ singularities. This 
leads to logarithmic divergency of the critical current 
( p^ ) in the limit T — > 0, similarly to the well known 
Riedel singularity of nonstationary supercurrent in SIS 
tunnel junctions at voltage eV = 2Ao, where the energy 
shift is due to the electric potential. At the same value 
of the exchange field h = the critical current 

changes its sign (i.e., the crossover from to tt contact 
occurs) for parallel magnetizations in the F layers [see 
Eq. (|l^)]. We emphasize that the scenario of the O-tt 
transition in our case differs from those studied before 
where the 7r-shift of the phase was either due to spatial 
oscillations of the order parameter in F layers or due to 
the proximity-induced phase rotation in S layers. In our 
case the phase does not change in either layer; instead, it 
jumps at the SF interfaces. This scenario is most clearly 
illustrated in the limit of large H where Eqs. (||), ( |l^ ) 
yield Fp oc — iAsgniJ whereas Fs oc A; thus the phase 
jumps by 7r/2 at each of the SF interfaces, providing 
the total TT-shift between Fpi{—H) and Fp2{H) [it is 
the phase difference between these two functions that 
determines the supercurrent according to Eq. (pT|)] . 

The considered effects take place only for sufficiently 
low I-barrier transparency. Indeed, it follows from Eq. 
( p^ ) that Gf{^) oc for small Q under condition 

h — If-fBH- As a result, the boundary condition (^) 
yields that at 



Q < min 



6 



7s 



dplBj' lB,I 



(16) 



the solutions ( |12D are not valid, since in this frequency 
range the effective transparency of the I interface (the 



parameter GfiGf2/ib,i [p^) increases and the spatial 
gradients in the F layers become large (the limit of large 
gradients is called "the KO-1 case" |lj, 0). In this case 
the nongradient term in Eq. (^) can be neglected and 
the general solution of the Usadel equation in the F lay- 



ers has the KO-1 form |13 : 

$ C - iM arctan [M {Bx + Q)] 



UJ 



1 



(17) 



where M = ^^(772 - 1) - C2, while C, B, Q and 
77 are integration constants. From Eqs. (|l|), ( [TtI ) 
it follows that the Green functions G, F and hence 
the contribution to the critical current from these fre- 
quencies are i?- independent. As a result, the bar- 
rier transparency parameter 7_b,7 provides the cutoff of 



(a) 



at 



the low-tcmpcrature logarithmic singularity of /, 
h = l/7sj\/ [see Eq. (0)]. According to Eq. (|l|), 
the critical current saturates at low temperature T* = 
Tcmin{^F /dpjB.i, 1b/ib.i)- We note that any asym- 
metry in the SFIFS junction will also lead to the cutoff 
of Ic"'^ divergency |l^ . The above estimates are done 
for the case of low barrier transparency, S,f / dp^Bj ^ 1 
and ^bIib.i <C 1. The opposite regime of high trans- 
parency deserves separate study. 

The general case. For arbitrary F-layer thick- 
nesses and interface parameters the boundary problem 
(0)"(|i) been solved numerically using iterative pro- 
cedure. Starting from trial values of the complex pair 
potentials A and the Green functions Gs,f we solve 
the resulting linear equations and boundary conditions 
for functions ^s,F- After that we recalculate Gs,f and 
A. Then we repeat the iterations until convergency is 
reached. The self-consistency of calculations is checked 
by the condition of conservation of the supercurrent ( p^ ) 
across the junction. We emphasize that our method is 
jully self-consistent: in particular, it includes the self- 
consistency over the superfluid velocity Vs, which is es- 
sential (contrary to the constriction case) in the quasi- 
one-dimensional geometry. The details of our numerical 
method will be presented elsewhere |]l9|| . 

Figure |l| shows Ic{F{) dependencies calculated at 
T = 0.05 Tc from the numerical solution of the bound- 
ary problem (|l|)-(^ for the fixed value of jbm = 1 
and a set of different F-layers thicknesses and the SF 
interface parameters 7. The normal junction resis- 
tance is Rn — Rb,i + 2i?B + 2pFdp/A. The curves 
dp/^F — are the limits of vanishing dp/^p ratio at 
fixed 7BAf and are calculated from Eqs. (p3|), (p^). For 
thin F layers the results depend only on the combina- 
tion 7j\,f — ^dp/^p. The enhancement of Ic and the 
crossover to the 7r-state are clearly seen for the antipar- 
allel and parallel orientations, respectively. In accor- 
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0.5 1.0 

Fig.|^. Enhancement of the critical current (antiparal- 
lel magnetizations, solid hnes) and the O-tt transition 
at which Ic changes its sign (parallel magnetizations, 
dashed lines) in the SFIFS junction at T/Tc = 0.05, 
7s A/ ~ 1, and jm = 0. Inset: the same for large values 
of 7JVI. 



1.5 



YsM= U y = 
r/r =0.02 




1.0 1.5 2.0 

Fig.^. Enhancement of the critical current (antiparal- 
lel magnetizations, solid lines) and the O-tt transition 
at which Ic changes its sign (parallel magnetizations, 
dashed lines) in the SFIFS junction: influence of tem- 
perature and barrier transparency. The dotted line cor- 
responds to T/Tc = 0.01 and S^f/cLf^bj ~ 0; the pa- 
rameters for other curves are given in the Figure. 



dance with the estimates given above, these effects take 
place for the values of the exchange field H close to 
■kTc- The enhancement disappears with increasing gra- 
dients in the F layers since the solution Eq. (|l^) loses 
its validity. This is illustrated in Fig.|l| by increasing 
the thickness dp or 7j\/. In particular, in the case of 
large 7m the enhancement is absent, in contrast to the 
statement in Ref. (see 111). 




Fig.^ Normalized density of states for spin "up" in the 
F layer for various exchange fields. 



Influence of temperature and barrier transparency 
on the critical current anomaly is shown in Figj^. One 
can see that, in accordance with the above estimate, 
the cutoff of /c°^ singularity is provided by finite tem- 
perature or barrier transparency. Namely, with the de- 
crease of the barrier strength parameter ^b.i the peak 
magnitude starts to drop when the ratio dp^B,i/^F be- 
comes comparable to T/Tc- With further decrease of 
dplB.i / S,F the singularity disappears, while the transi- 
tion to the TT-state shifts to large values of H . 

Figure ^ demonstrates the DoS in the F layers for 
one spin projection, calculated numerically in the limit 
of small I-barrier transparency. At i7 = we reproduce 
the well-known minigap existing in SN bilayer. At fi- 
nite H the gap shifts in energy (asymmetrically) and 
the peak in the DoS reaches zero energy at h — 1/jbm- 
One can see that even for a small value 7m = 0.05 the 
peaks are rather broad, this is the reason why the sin- 
gularity in /i"'' is suppressed by 7m very rapidly. 

In the practically interesting limit of finite F-layer 
thickness (see Fig.^ the numerical calculations show 
monotonic suppression of Ic with increase of the ex- 
change field H for antiparallel magnetizations of the F 
layers and the O-tt crossover for the parallel case. One 
can see from Fig.^ that for given temperature and thick- 
ness of the F layers it is possible to find the value of the 
exchange field at which switching between parallel and 
antiparallel orientations will lead to switching of Ic from 
nearly zero to a finite value (or to switching between 
and TT states). This effect may be used for engineer- 
ing cryoelectronic devices manipulating spin-polarized 
electrons. 
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Fig.y. Critical current in SF1F2S junction: switching 
effect. T/Tc — 0.5, the solid and dashed lines corre- 
spond to the antiparallel and parallel orientations of 
magnetizations, respectively. Inset: thermally induced 
0-7r crossover in the parallel case. 

The case of parallel F-layers magnetizations in the 
absence of the I barrier corresponds to the standard SFS 
junction where the O-tt transition is possible due to spa- 
tial oscillations of induced superconducting ordering in 
the F layer. The thermally induced O-tt crossover in SFS 
junction was observed in Ref. where simple theory 
based on the linearized Usadel equations was also pre- 
sented. Here we show such a crossover (see the inset in 
Fig.^ from the fully self-consistent solution in the range 
of the exchange fields corresponding to that of Ref. . 
Comparison to the experimental data and more detailed 
results of our model will be given elsewhere (l9) . 

In conclusion, we have presented a general method 
to solve the Usadel equations in SFIFS jimctions self- 
consistently. Using our method, we have investigated 
theoretically the Josephson current in SFIFS and SFS 
junctions as a function of relative F-layers magnetiza- 
tions, thicknesses and parameters of the S/F and F/F 
interfaces. We have identified the physical mechanisms 
of the critical current enhancement and of the O-tt tran- 
sition in these junctions. 
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